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Abstract 

In the framework of a variational method with a single trial function an accurate study of the 
lowest gerade l g and ungerade l u electronic states of the molecular ion i7^~ in a magnetic field 
is performed. Magnetic field ranges from to 4.414 x 10 13 G and orientations of the molecular 
axis with respect to the magnetic line 0° < 6 < 90° are considered. A one-parameter gauge 
dependent vector potential is used in the Hamiltonian, which is finally variationally optimized. 
A well pronounced minimum on the total energy surface of the (ppe) system in both l g and l u 
states is found for all magnetic fields and orientations studied. It is shown that for both states the 
parallel configuration (6 = 0) at equilibrium always corresponds to the minimal total energy. It 
is found that for a given 9 for both states the magnetic field growth is always accompanied by an 
increase in the total and binding energies as well as a shrinking of the equilibrium distance. We 
demonstrate that for B > 1.8 X 10 11 G the molecular ion can dissociate, H^ — ► H + p, over a 
certain range of orientations (6 cr < 6 < 90°), where the minimal 8 cr ~ 25° occurs for the strongest 
magnetic field studied, B = 4.414 x 10 13 G. For B < 10 12 G the ion H^ in l g , l u states is the most 
compact, being in the perpendicular configuration (0=90°), whereas for B > 10 12 this occurs for 
an angle < 90°. For the l g state in any orientation, with the magnetic field growth at B ~ 10 11 
G, a two-peak electronic distribution changes to single-peak one. 
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I. INTRODUCTION 

There are theoretical qualitative indications that in the presence of a strong magnetic field 
the physics of atoms and molecules gets a wealth of new, unexpected phenomena even for the 
simplest systems Jl], @. In practice, the atmosphere of neutron stars, which is characterized 
by the existence of enormous magnetic fields 10 12 — 10 13 G, provides a valuable paradigm 
where this physics could be realized. 

One of the first general features observed in standard atomic and molecular systems 
placed in a strong magnetic field is an increase in the binding energy, accompanied by a 
drastic shrinking of the electron localization length. It leads to a decrease of the equilibrium 
distance with magnetic field growth. This behavior can be considered as a consequence of the 
fact that at large magnetic fields the electron cloud takes a needle-like form extended along 
the magnetic field direction and the system becomes effectively quasi-one-dimensional. Fur- 
thermore the phenomenon of quasi-one-dimensionality enhances the stability of standard 
atomic and molecular systems. They become elongated along the magnetic line forming 
molecules of the type of linear polymers (for details see review papers || £|). It also hints 
at the existence of exotic atomic and molecular systems which do not exist in the absence 
of a magnetic field. Motivated by these simple observations it was shown in || |(| that 
exotic one-electron molecular systems H% + and H^ ++ can exist in sufficiently strong mag- 
netic fields in the form of linear polymers. However, the situation becomes much less clear 
(and also much less investigated) when the nuclei are not aligned with the magnetic field 
direction, thus not forming in general a linear system. Obviously, such a study would be 
important for understanding the kinetics of a gas of molecules in the presence of a strong 
magnetic field. As a first step towards such a study, even the simplest molecules in dif- 
ferent spatial configurations deserve attention. The goal of the present work is to attempt 
to make an extensive quantitative investigation of the lowest electronic states of H^~ in a 
(near) equilibrium position in the framework of a single approach in its entire complexity: 
wide range of magnetic field strengths (0 — 4 x 10 13 G) and arbitrary orientation of the 
molecular axis with respect to magnetic line. 

It is well known that the molecular ion H^~ is the most stable one-electron molecular 



system in the absence of a magnetic field. It remains so in the presence of a constant magnetic 

l 3 



field as well, but it turns out that for a magnetic field B > 10 13 G the exotic ion Hjj" + appears 



to be the most bound (see 0). The ion H^" has been widely studied, both with and without 
the presence of a magnetic field, due to its importance in astrophysics, atomic and molecular 
physics, solid state and plasma physics (see P-|22f and references therein). The majority of 
the previously performed studies was focused on the case of the parallel configuration, when 
the angle between the molecular axis and the magnetic field direction was zero, 8 = 0, with 



an exception of [14" ], where a detailed quantitative analysis was performed for any 9 but 
for B = la.u. Previous studies were based on various numerical techniques, most of all - 
different versions of the variational method, including the Thomas-Fermi approach. As a 
rule, in these studies the nuclear motion was separated from the electronic motion using 
the Born-Oppenheimer approximation at zero order - assuming protons as infinitely heavy 
charged centers. It was observed at a quantitative level that magnetic field growth is always 
accompanied by an increase in the total and binding energies, as well as shrinking of the 
equilibrium distance. As a consequence it led to a striking conclusion about the drastic 
increase in the probability of nuclear fusion for H^~ in the presence of a strong magnetic field 

In the present study we perform accurate calculations for the lowest l g and l u electronic 
states of Hj in (near)equilibrium position for magnetic fields B = — 4.414 x 10 13 G and 
arbitrary orientation of the molecular axis towards the magnetic line. Since the Hamiltonian 
is gauge-dependent a choice of the form of the vector potential is one of the crucial points 
in our study. We construct state-of-the-art, 'adequate', trial functions with a variationally 
optimized gauge dependence consistent with the choice of vector potential. Although an ap- 
propriate position of the gauge origin|27| may be important (especially for large internuclear 



distances, this lies beyond the scope of the present article and will be discussed elsewhere) 
we place the gauge origin in the middle between the two nuclei (charged centers) and keep 
it fixed. For the parallel configuration the present work can be considered as an extension 
(and also an improvement) of previous work | 20fl . It is necessary to emphasize that we en- 



countered several new physical phenomena which occur when the molecular axis deviates 
from the magnetic field direction. In particular, if the magnetic field is sufficiently strong, 
B ~ 10 11 G the ion H^ can dissociate to H+p [pS|j . This means that even though the positive 
binding energy of Hj in the optimal configuration ensures its existence, even under the high 
temperature conditions prevailing on the surface of neutron stars (10 — 100 eV), there is a 
certain probability of dissociation. The behavior of the equilibrium distance as a function 



of 9 reveals another surprising feature: for magnetic fields B< 10 12 G the molecule is most 
compact in the perpendicular configuration, while for B > 10 12 G this occurs for a certain 
angle 9 < 90°. We find that for the l g state at any orientation in the weak field regime the 
electronic distribution peaks at the positions of the protons, while at large magnetic fields 
the electronic distribution is characterized by one maximum which occurs at the midpoint 
between two protons. This change appears around B ~ 10 10-11 G with a slight dependence 
on the inclination angle 9. From a physical point of view the former means that the electron 
prefers to stay in the vicinity of a proton. It can be interpreted as a dominance of the 
if -atom plus proton interaction. The latter situation implies that the electron is 'shared' by 
both protons and hence a separation to H-aXom. plus proton cannot be done. Therefore, we 
can call the two-peak situation an ionic coupling, while the one-peak case is assigned to the 
covalent coupling, although this definition differs from the one widely accepted in textbooks 
(see for example @]). Thus, we can conclude that a new phenomenon appears - as the 
magnetic field grows the type of coupling changes from ionic to covalent. 

Atomic units are used throughout (h=m e =e=l) albeit energies are expressed in Rydbergs 
(Ry). In particular, the magnetic field B is given in a.u. with Bq = 2.35 x 10 9 G. 

II. THEORY 

The Hamiltonian which describes the Hjj~ molecular ion placed in a uniform constant 



magnetic field directed along the z-axis, B = (0, 0, B) is given by (see, for example, p3 



H = f + ?------(pA)+A\ (1) 

R ri r 2 

(see Fig.l for notations), where p = —iV is the electron momentum, A is a vector potential 
corresponding to the magnetic field B. 

The vector potential is gauge-dependent and hence it is defined ambiguously. Therefore 
the Hamiltonian (JI|) is also explicitly gauge-dependent although the energies and other ob- 
servables do not. Thus, when solving the eigenvalue problem corresponding to the above 
Hamiltonian within an approximate method, the energies (as well as other observables) will 
in general be gauge- dependent (only the exact ones would be gauge-independent). Hence, 
one can choose the form of the vector potential in a certain optimal way, for instance, which 
would lead to minimal approximate total energy of the ground state. In order to realize this 




FIG. 1: Geometrical setting for the H^ ion in a magnetic field directed along the z-axis. 

idea let us take a certain one-parametric family of vector potentials corresponding to the 
constant magnetic field B (e.g. see M) 



A = B((Z-l)y,£x,0) 



(2) 



where £ is the parameter to be chosen in an optimal way. The well known symmetric (or 
circular) and asymmetric gauges are particular cases of (0) corresponding to £=1/2 and 
£ = 1, respectively (see, for instance, [p3l ). By substituting (H) into (|]) we arrive at the 
Hamiltonian 
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H = -V 2 + ----- + iB[(£ - l)yd x + ixd y \ + B% l x l + (1 - o 2 y 



(3) 



The idea of choosing an optimal gauge is widely exploited in quantum field theory. Usu- 
ally, the gauge is fixed following a certain convenience criterion or for technical simplicity. 
Regarding our present problem of H^~ in a magnetic field, there were also some attempts 
to discuss the gauge fixing (see, for instance, [[H] and references therein). Perhaps the 
first constructive (and remarkable) attempt was realized by Larsen || in his study of the 
ground state. In particular, he explicitly showed that gauge dependence of the variational 
results can be quite significant and even an oversimplified optimization procedure improves 
drastically the accuracy of the numerical results. 

It is rather evident that there exists a certain gauge for which the Hamiltonian (|]) has the 
ground state eigenfunction given by a real function [pSfl . Thus, we can use real trial functions 
with explicit dependence on the gauge parameter £. This parameter is fixed by performing 
a variational optimization of the energy. Therefore, as a result of the minimization we find 



both a variational energy and a gauge for which the eigenfunction is real, as well as the 
corresponding Hamiltonian. One can show that for a system possessing axial (rotational) 
symmetry |3D] the optimal gauge is the symmetric one £ = 1/2. It is precisely this gauge 



which was overwhelmingly used (without any explanations) in the majority of previously 
performed H;j" studies in the parallel configuration [3-22,25]. However, this is not the case 
if 9 7^ 0°. For the symmetric gauge the exact eigenfunction now becomes complex, hence 
complex trial functions should be used. But we can avoid complexity by adjusting the gauge 
so as to continue with real eigenfunctions. This justifies the use of real trial functions. Our 
results (see Section 3) lead to the conclusion that for both studied states the optimal gauge 
parameter is £ G [0.5, 1]. In the limit B 3> leu., the parameter £ tends to one corresponding 
to asymmetric gauge. Likely, this tendency will continue for other excited states. 

One can easily see that the expectation value of the term ~ B in Eq. (Q) vanishes when it is 
calculated with respect to any real normalizable function. Hence, without loss of generality, 
this term can be omitted in the Hamiltonian (fj). It gives an essential simplification of the 
computational procedure. Our choice of variational trial function \l/o can be formulated as 
follows J22|: we construct a real trial function for which the potential Vq = -^ reproduces 



the original potential near Coulomb singularities and its growing terms at large distances; it 
also supports all symmetries of the problem at hand and the gauge parameter £ is included 
into the trial function explicitly. This recipe (for symmetric gauge) was successfully exploited 
in the previous study of H^~ in strong magnetic fields in the parallel configuration |2U| . This 
recipe was also used to construct trial functions when the question of the existence of the 
exotic ions Hj~ + and H^ ++ is considered |)J [| . 

The Hamiltonian ([!]) gives rise to different symmetry properties depending on the orienta- 
tion of the magnetic field with respect to the internuclear axis. The most symmetric situation 
corresponds to 9 = 0°, where invariance under permutation of the (identical) charged centers 
P : (1 <-> 2) together with P z : (z — > —z) symmetry hold. Since the angular momentum 
projection C. z = m is conserved, P z accounts also for the degeneracy m — > —m. Thus, we 
classify the states as lcr gtU , 2a 9iU , . . . l7T gtU , 27r S)U . . . lS 9iU , 25 9tU . . . , where the numbers 1,2,... 
refer to the electronic states in increasing order of energy. The labels a, n, 5 . . . are used 
to denote \m\ = 0,1,2..., respectively, the label g (u) gerade (ungerade) is assigned to 
the states of even (odd) parity P of the system. At 9 = 90° the Hamiltonian still remains 
invariant under the parity operations P and P z , while the angular momentum projection 



is no longer conserved and m is no longer a quantum number. The classification in this 
case is lf u , 2^ u , . . . , where sign +(— ) is reserved to denote even (odd) z-parity. Eventually, 
for arbitrary orientation the only parity under permutations P is conserved. In general, we 
refer to the lowest gerade and ungerade states in our study as l g and l u . This is the only 
notation which make sense for all orientations 0° < 9 < 90°. 

One of the simplest trial functions for l g state which meets the requirements of our 
criterion of adequacy is 

^ = e -^(r 1 +r 2 ) e -B[/3 lx ^x 2 +/3 1 y(l-Oy 2 ] (4) 



(cf. p0| , |22f), where «i, (3\ x and (3\ y are variational parameters and £ is the parameter 
of the gauge (2). The first factor in this function, being symmetric under interchange of 
the charged centers r± <-» r 2 , and corresponding to the product of two ls-Coulomb orbitals 
centered on each proton, is nothing but the celebrated Heitler-London approximation for the 
ground state la g . The second one is the lowest Landau orbital corresponding to the vector 
potential of the form Eq. (^|). So, the function (|j) can be considered as a modification of the 
free field Heitler-London function. Following the experience gained in studies of H^ without 
magnetic field it is natural to assume that Eq. (^) is adequate to describe internuclear 
distances near equilibrium. This assumption will be checked (and eventually confirmed) a 
posteriori, after making concrete calculations (see Section 3). 

The function (H) corresponds to an exact eigenfunction of the potential 
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The potential V\ reproduces the functional behaviour of the original potential (|3|) near 
Coulombic singularities and at large distances. These singularities are reproduced exactly 
when (3\ x = f3\ y = 1/2 and at\ = 1. 
Another trial function meets the requirements of our criterion of adequacy as well 

v|> 2 = (e~ a2n + Pe~ a ' 2T2 \ e - B \P**t a ?+fov(. 1 -0v 2 ] ; (5) 



(cf. [|0|, |22fl). It is the celebrated Hund-Mulliken function of the free field case multiplied by 
the lowest Landau orbital, where a 2 , fax and fa y are variational parameters. The parameter 
P has the meaning of the parity, being equal to +1 for l g and —1 for l n . From a physical 
point of view this function describes the interaction between a hydrogen atom and a proton 
(charged center), and also models the possible decay mode of H^" into hydrogen atom plus 
proton. Thus, one can naturally expect that for sufficiently large internuclear distances R 
this function prevails giving a dominant contribution. Again this assumption will be checked 
a posteriori, by concrete calculations (see Section 3). 

There are two natural ways to incorporate the behavior of the system both near equilib- 
rium and at large distances into a single trial function. It is to make a linear or nonlinear 
interpolation. The linear interpolation is given by a linear superposition 

q 3a = A 1 q 1 + A 2 ty 2} (6) 

where Ai or A 2 are parameters and one of them is kept fixed by a normalization condition. 
The function @ makes sense for l g state only since the function (||) does not exist for l u 
state. Therefore, conditionally, for state l u we simply put the parameter Ai = 0. In turn, 
the simplest nonlinear interpolation is of the form 

\£n h = ( e - a 3»'i-Q4r'2 _|_ p e -a 3 r 2 -a 4 r 1 \ e -B[p 3x £x 2 +p 3y (l-£)y 2 ] /y\ 

(cf. [F2U], |22[]), where a 3 , a 4 , fa x and fa y are variational parameters and P = ±1 is the parity 
of the state. This is a Guillemin-Zener function of the free field case multiplied by the lowest 
Landau orbital. If a 3 = a 4 , the function (|7|) coincides with (|J). If a 4 = 0, the function (^) 
coincides with (||). 

The most general Ansatz is a linear superposition of trial functions @ and (|7|), 

^ = Ai^i + A 2 ^ 2 + A 3 ^ 3b , (8) 

where we fix one of the A's and let all other parameters vary. For the l u state we put the 
parameter Ai = since the function (|J) does not exist for this state. Finally, the total 
number of variational parameters in @, including R and £, is fourteen for the l g state 
and ten for the l u state, respectively. For the parallel configuration, 9 = 0°, the parameter 
£ = 0.5 and hence fa x = fa yi fa x = fay, fax = fay, an d the number of free parameters is 
reduced to ten for the l g state and seven for the l u state, respectively. Finally, with the 



function (§) we intend to describe the lowest (un)gerade state for all magnetic fields where 
non-relativistic consideration is valid for B < 4.414 x 10 13 G, and for all orientations of the 
molecular axes. 

Our variational calculations were performed by using the minimization package MINUIT 
of CERN-LIB. Numerical integrations were carried out with relative accuracy ~ 10~ 7 by use 
of adaptive NAG-LIB (D01FCF) routine. All calculations were performed on a Pentium-Ill 
PC 750MHz. 

III. RESULTS AND DISCUSSION 

In this Section we present the results for the l g and the l u states of Hj for magnetic fields 
ranging from B = G through B = 4.414 x 10 13 G, where a non-relativistic consideration 
is relevant and relativistic corrections can be neglected (see Jl9[ for discussion), and for 
orientations ranging from 0° (parallel configuration) to 90° (perpendicular configuration). 

A. l g state 

Before coming to a concrete quantitative consideration we must state that for all explored 
magnetic fields, B = — 4.414 x 10 13 G, and all orientations we found that a well pronounced 
minimum in the total energy of the system (ppe) appears at finite internuclear distance. 



This minimum always corresponds to positive binding energy Ef, [31] and hence it implies 
the existence of the H^" molecular ion for magnetic fields B = — 4.414 x 10 13 G and any 
orientation of the molecular axis. This is at variance with the statement by Khersonskij 
J7J about the non-existence of a minimum at finite distances on the total energy surfaces 
at sufficiently strong magnetic fields for some, far from parallel orientations and therefore 
non-existence of the molecular ion H^~. Presumably, this statement of non-existence is an 
artifact arising from an inappropriate choice of trial functions and therefore a consequent loss 
of accuracy. It is worth emphasizing that the variational study in was carried out with 
a trial function somewhat similar to that of Eq.@, but which does not fulfill our criterion 
of adequacy in full. The potential corresponding to this function reproduces correctly the 
original potential near Coulomb singularities and ~ p 2 -growth at large distances. However, 
it generates growing terms ~ p which implies a reduction in the rate of convergence of 
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a perturbation theory for which the variational energy represents the first two terms (see 
discussion in [f24|). 

In Tables I, II and III the results for the total energy (Et), binding energy (E^) and 
equilibrium distance {R eq ) are shown for 6=0°, 45° and 90°, respectively. As seen in Table 
I, our results for 9 = 0° lead to the largest binding energies for B > 10 11 in comparison with 
other calculations. As for B < 10 11 G, our binding energies for the parallel configuration 
appear to be very close (of the order of < 10~ 4-5 in relative deviation) to the variational 



results at Wille [[1211 , which are the most accurate so far in this region of magnetic field 
strengths. They are based on a trial function in the form of a linear superposition of < 500 
Hylleraas type functions. It is quite amazing that our simple trial function (8) with ten 
variational parameters gives comparable (for B < 10 11 G) or even better (for B > 10 11 
G) accuracy. It is important to state the reason why the trial function Jl2[ fails being, 
increasingly inaccurate as a function of magnetic field growth for B > 10 11 G. An explanation 
of this inaccuracy is related to the fact that in the (x, y)- directions the exact wave function 
decays asymptotically as a Gaussian function, unlike Hylleraas functions which decay like the 
exponential of a linear function. The potential corresponding to the function EJ reproduces 
correctly the original potential near Coulomb singularities but fails to reproduce ~ p 2 -growth 
at large distances. It implies a zero radius of convergence of the perturbation theory for 
which the variational energy represents the first two terms (see discussion in [|24}|). 

A failure to adequately reproduce the asymptotic behavior of the wave functions leads 
to a failure in the proper description of the electronic cloud shrinking in transversal to 
magnetic line directions which become more and more essential with magnetic field growth. 
The origin of this shrinking is the Lorentz force action. The above drawback can easily be 
fixed by modifying the Hylleraas function by multiplication on the lowest Landau orbital. 
Nevertheless, it is worth emphasizing that it is quite surprising that a linear superposition 
of Hylleraas type functions, which form a natural basis for the free-field case, still allows one 
to get high accuracy results for magnetic fields B < 10 11 G. Implicitly, it indicates that the 
ground state wave function for these magnetic fields does not deviate drastically from free- 
field behavior. It is obvious that taking the above-mentioned modified Hylleraas functions 



in the analysis made by Wille [F3 would allow us to reach the same accuracies but with 
a much smaller basis. We can treat the results by Larsen || as an explicit demonstration 
that an insertion of the Landau orbitals (multiplied by a Gaussian in z function) in the 
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TABLE I: 


Total Er, 


binding Ej, energies and equilibrium dist 


ance R eq 


for the state l g in parallel 


configuration, (9 = 0° 


t This value is 


taken from [M 






13 




E T (Ry) 


E b (Ry) 


R eq (a.u.) 




B = 




-1.20525 


— 


1.9971 


Present^ 






-1.20527 


— 


1.997 


Willc H 


10 9 G 




-1.15070 


1.57623 


1.924 


Present 






-1.15072 


1.57625 


1.924 


Willc [[| 


1 a.u. 




-0.94991 


1.94991 


1.752 


Present 






— 


1.9498 


1.752 


Larsen [EI 






-0.94642 


1.94642 


1.76 


Kappes et al Ej 


10 10 G 




1.09044 


3.16488 


1.246 


Present 






1.09031 


3.16502 


1.246 


Willc [Q 


10 a.u 




5.65024 


4.34976 


0.957 


Present 






— 


4.35 


0.950 


wilic Hi 






— 


4.35 


0.958 


Larsen [EI 






— 


4.3346 


0.950 


Vincke et al (lo) 


10 11 G 




35.0434 


7.50975 


0.593 


Present 






35.0428 


7.5104 


0.593 


Willc H 






— 


7.34559 


0.61 


Lai et al Jli| 


100 a.u. 




89.7096 


10.2904 


0.448 


Present 






— 


10.2892 


0.446 


Wille H 






— 


10.1577 


0.455 


Wunner et al |21[ 






— 


10.270 


0.448 


Larsen [El 






— 


10.2778 


0.446 


Vincke et al Jlo| 


10 12 G 




408.3894 


17.1425 


0.283 


Present 






— 


17.0588 


0.28 


Lai et al |l8) 






408.566 


16.966 


0.278 


Willc [0 


1000 a.u 




977.2219 


22.7781 


0.220 


Present 






— 


21.6688 


0.219 


wilic Hi 






— 


22.7069 


0.221 


Wunner et al |jlt 






— 


22.67 


0.222 


Larsen [EI 






— 


22.7694 


0.219 


Vincke et al Jltj 


10 13 G 




4219.565 


35.7539 


0.147 


Present 






4231.82 


23.52 


0.125 


wilic h 






— 


35.74 


0.15 


Lai et al Jli| 


4.414 x 10 13 


G 


18728.48 


54.4992 


0.101 


Present 



trial function (hence not in a fully adequate way as follows from our recipe) together with 
an optimization of the gauge dependence (in a different manner than what we propose) 
leads to rather accurate results. Finally, it is worth mentioning that in the domain of very 
strong magnetic fields, B > 10 12 G, our results are more accurate than those obtained by 



the Thomas-Fermi method 18 . 
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TABLE II: Total Et, binding Ef, energies and equilibrium distance R eq for the l g state at 9 = 45°. 
Optimal value of the gauge parameter £ is given (see text). 



B 


E T (Ry) 


E b (Ry) 


Req (a.U.) 


£ 


10 9 G 


-1.14248 


1.56801 


1.891 


0.5806 


1 a.u. 


-0.918334 


1.91833 


1.667 


0.5846 


10 10 G 


1.26195 


2.99337 


1.103 


0.5958 


10 a.u. 


6.02330 


3.97670 


0.812 


0.6044 


10 11 G 


36.15638 


6.39681 


0.466 


0.6249 


100 a.u. 


91.70480 


8.29520 


0.337 


0.6424 


10 12 G 


413.2988 


12.2331 


0.198 


0.6894 


1000 a.u. 


985.1956 


14.8044 


0.147 


0.7151 


10 13 G 


4236.342 


18.9769 


0.097 


0.8277 


4.414 x 10 13 G 


18760.77 


22.2047 


0.073 


0.9276 



The results for 9 = 45° are shown in Table II, where a gradual shortening of the equi- 
librium distance is accompanied by an increase of total and binding energies with magnetic 
field growth. It is worth noting that the parameter £ evolves from about 0.5 to 0.93 with 
magnetic field growth, thus changing from symmetric gauge for weak fields to an almost 
asymmetric one for strong ones. This phenomenon takes place for all orientations 8^0, 
becoming more and more pronounced with the inclination angle growth (see below). We are 
unaware of any other calculations for 9 = 45° to compare with. 

For the perpendicular configuration (9 = 90°) the results are presented in Table III. 
Similar to what appeared for the parallel configuration (see above) our results are again 
slightly less accurate than those of Wille for B < 10 10 G becoming the most accurate for 
stronger fields. In particular, it indicates that a domain of applicability of the trial function, 
taken in the form of a superposition of Hylleraas type functions, reduces when the inclination 



grows. The results reported by Larsen || and by Kappes-Schmelcher [14] are slightly worse 



than ours although the difference is very small. The evolution of the gauge parameter follows 
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TABLE III: Total E T , binding E b energies 
perpendicular configuration, 6 = 90°. Optimal 



distance R eq for the l g state in 



and equilibrium 

value of the gauge parameter £ is given (see text). 



B 


E T (Ry) 


E b (Ry) 


Req (a.U.) 


e 




10 9 G 


-1.137342 


1.56287 


1.875 


0.6380 


Present 






1.56384 


1.879 




Wille || 


1 a.u. 


-0.89911 


1.89911 


1.636 


0.6448 


Present 




— 


1.8988 


1.634 




Larsen || 




-0.89774 


1.8977 


1.65 




Kappes et al 14 1 


10 10 G 


1.36207 


2.89324 


1.059 


0.6621 


Present 




— 


2.8992 


1.067 




Wille [12] 


10 a.u. 


6.23170 


3.76830 


0.772 


0.6752 


Present 




— 


3.7620 


0.772 




Larsen || 


10 11 G 


36.7687 


5.78445 


0.442 


0.7064 


Present 




— 


5.6818 


0.428 




Wille [12] 


100 a.u. 


92.7346 


7.26543 


0.320 


0.7329 


Present 




— 


7.229 


0.320 




Larsen || 


10 12 G 


415.5621 


9.96986 


0.196 


0.8034 


Present 




— 


4.558 


0.148 




Wille 0| 


1000 a.u. 


988.3082 


11.6918 


0.151 


0.8520 


Present 




— 


11.58 


0.1578 




Larsen |J 


10 13 G 


4241.470 


13.8490 


0.113 


0.9359 


Present 


4.414 x 10 13 G 


18767.50 


15.4700 


0.0937 


0.9795 


Present 



a similar trend as was observed at 9 — 45°, varying from £ = 0.64 to £ = 0.98 with magnetic 
field growth from B = 10 9 G to B = 4.414 x 10 13 G [g| . 



The total energy dependence of H^~ as a function of the inclination angle 9 for different 
magnetic fields is shown in Fig.2. The dotted line corresponds to the if -atom total energy 
in the magnetic field. For weak magnetic fields the hydrogen atom total energy is higher 
than that of the H^-ion one. However, for B > 1.8 x 10 11 G the total energy of the if -atom 
becomes lower than the total energy of the H^~-ion for some angles 9 cr < 9 < 90°. It implies 
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that one proton in the H^ system can go to infinity. Thus, it manifests the appearance of 
a dissociation channel H^ — > H + p. At first, the dissociation occurs at 9 = 90°. Afterward 
the domain of inclinations with allowed dissociation widens with the magnetic field growth 
reaching 25° < 9 ^ 90° for B = 4.414 x 10 13 G. The dependence of the critical angle 9 cr on 
the magnetic field is shown in Fig.3. Naively, it looks like the rate of dissociation is maximal 
at 9 = 90°. However, a precise conclusion depends on the form of the barrier or, in other 
words, on the form of the potential surface in 9, R. The rate of dissociation as a function 
of inclination angle and magnetic field is not studied in detail in present work and will be 
published elsewhere. 

We observe that for fixed value of the magnetic field strength B the binding energy Ej, 
as a function of 9 always decreases when changing from the parallel to the perpendicular 
configuration (see Fig. 2). A similar picture holds for all studied values of the magnetic field 
strength. Thus, we can draw the conclusion that the molecular ion becomes less and less 
stable monotonically as a function of inclination angle growth. It confirms the statement 
|7], |9|, 0, PH , that the highest molecular stability of the l g state occurs for the parallel 
configuration. We extend its validity to magnetic field strengths B < 4.414 x 10 13 G. It is 
worth emphasizing that the rate of increase of binding energy with magnetic field growth 
depends on the inclination - it slows down with inclination increase. This effect means that 
H2 in the parallel configuration becomes more and more stable towards rotations - the 
energy of the lowest rotational state should increase rapidly with magnetic field (see Table 
V and discussion there). 

Regarding the internuclear equilibrium distance R eq , one would straightforwardly expect 
that it will always decrease with inclination growth. Indeed, for all studied magnetic fields 
we observe that R eq at 9 = 0° is larger than R eq at 9 = 90° (cf. Tables I, III). This can be 
explained as a consequence of the much more drastic shrinking of the electronic cloud in the 
direction transverse to the magnetic field than in the longitudinal one. Actually, for magnetic 
fields B < 10 12 G the equilibrium distance R eq decreases monotonically with inclination 
growth, as seen in Fig.4. However, this trend breaks down for higher magnetic fields where 
the shortest equilibrium distances occur for orientations 9 m i n < 90° (!). Furthermore, as the 
magnetic field grows the molecular ion becomes the most compact for smaller angles, being 
Omin ~ 60° for B < 10 12 G and then going down to 9 min ~ 30° for B = 4.414 x 10 13 G. 
The minimal value of R eq also deepens in comparison with R eq at 9 = 90°. For example, 
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FIG. 2: H2 total energy (Et) for the ground state l g as function of the inclination angle 8 for 



different magnetic fields. The dotted lines correspond to the H-atom total energy taken from [18]. 



for B = 4.414 x 10 G it becomes almost twice as small as R Pn at 9 = 90° . All these 



>'<1 



irregularities appear when at the same time the binding energy decreases monotonically as 
inclination grows (see Fig. 2). We do not have a physical explanation of this phenomenon 
yet. 
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FIG. 3: HZ : domains of dissociation <-» non-dissociation for the l n state. 



Fig. 5 illustrates the above-mentioned non-monotonic behaviour of the internuclear dis- 
tances near the equilibrium position, R < 0.3 a.u., for different orientations at B — 10 13 G. 
For all values of 9 a clear minimum in R develops. For certain orientations (8 > 45°) the 
potential energy curves are situated above the if -atom total energy. The total energy curves 
lie at increasingly higher energies as inclination grows. Hence, one can draw the conclusion 
that the growth of the inclination angle leads to an excitation of the H^ system. It is quite 
interesting to make a comparison of the present situation with what appears in chemistry. 
The typical situation for the behaviour of molecular electronic excitations of standard chemi- 
cal systems (in the absence of a magnetic field) is characterized by an increase of equilibrium 
distance with energy growth. This is explained either by enhancement of the antibonding 
character or suppression of the bonding character. Our situation is the opposite - in general 
the antibonding character is suppressed or, equivalently, the bonding character is enhanced. 
However, for strong magnetic fields B > 10 12 G, where an abnormal behavior of the equilib- 
rium distance is seen (see Fig. 4), the situation is reminiscent of the chemical one. Starting 
from 9 cr the growth of inclination leads to an increase in the equilibrium distance. We do 
not know if this analogy makes physical sense or it is simply a coincidence. 

In order to characterize the electronic distribution for different orientations we have cal- 
culated the expectation values of the transversal < p > and longitudinal < \z\ > sizes of the 
electronic cloud (see Table IV). The ratio 

<P> 



< \z\ > 



< 1 
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1.925 
1.915 

1.905 
R sq (a.u.) 
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FIG. 4: H.2 equilibrium distance as function of the inclination angle 9 for the l g state. 

quickly decreases with magnetic field growth, especially for small inclination angles. It 
reflects the fact that the electronic cloud has a more and more pronounced needle-like form, 
oriented along the magnetic line as was predicted in [T|-(2|]. The behaviour of < p > itself 
does not display any unusual properties, smoothly decreasing with magnetic field, quickly 
approaching the cyclotron radius for small inclinations at large magnetic fields. On the 
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4255 




4215 



0.1 0.2 

R (a.u.) 



0.3 



FIG. 5: H^ potential curves (total energy Ex vs R) for the l g state at B = 10 13 G and different 
inclinations. The dotted line corresponds to the H-atom total energy taken from fll8|| . Positions of 
minimum are marked by bullets. 

contrary, the < \z\ > behaviour reveals some surprising features. At the beginning, at small 
magnetic fields, the < \z\ > expectation value monotonically decreases with inclination, but 
then after some irregular behavior at 10 12 < B < 10 13 , it begins a monotonic increase with 
inclination. It is quite striking that there is a domain of magnetic fields where < \z\ > has 
almost no dependence on inclination. 

As a result of our analysis the parallel configuration turned out to be optimal for all stud- 
ied magnetic fields. Therefore, it makes sense to perform a study of the lowest vibrational 
state and also the lowest rotational state (see Table V). In order to do this we separate 
the nuclear motion along the molecular axis near equilibrium in the parallel configuration 
(vibrational motion) and deviation in 9 of the molecular axis from 6 = 0° (rotational mo- 
tion). The vicinity of the minimum of the potential surface E(8,R) at 9 = 0°,R = R eq 
is approximated by a quadratic potential and hence we arrive at a two-dimensional har- 
monic oscillator problem in (R,6). Corresponding curvatures near the minimum define the 
vibrational and rotational energies (for precise definitions and discussion see, for example, 
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TABLE IV: l g state: Expectation values of the transversal < p > and longitudinal < \z\ > sizes 
of the electron distribution in H^ in a.u. at different orientations and magnetic field strengths. At 
6 = 0° the expectation value < p > almost coincides to the cyclotron radius of electron. 



B 




<p> 






< \ z \ > 






0° 


45° 


90° 


0° 


45° 


90° 


10 9 G 


0.909 


1.002 


1.084 


1.666 


1.440 


1.180 


1 a.u. 


0.801 


0.866 


0.929 


1.534 


1.313 


1.090 


10 10 G 


0.511 


0.538 


0.569 


1.144 


0.972 


0.848 


10 a.u. 


0.359 


0.375 


0.396 


0.918 


0.787 


0.708 


10 11 G 


0.185 


0.193 


0.205 


0.624 


0.542 


0.514 


100 a.u. 


0.123 


0.129 


0.139 


0.499 


0.443 


0.431 


10 12 G 


0.060 


0.065 


0.074 


0.351 


0.324 


0.340 


1000 a.u. 


0.039 


0.043 


0.054 


0.289 


0.275 


0.290 


10 13 G 


0.019 


0.025 


0.037 


0.215 


0.221 


0.256 


4.414 x 10 13 G 


0.009 


0.017 


0.030 


0.164 


0.191 


0.232 



H). We did not carry out a detailed numerical analysis, making rough estimates of the 
order of 20%. For example, at B= 10 12 G we obtain E v n, = 0.276 Ry in comparison with 
E V ib = 0.259 Ry given in ||, where a detailed variational analysis of the potential electronic 
curves was performed. Our estimates for the energy, E vi b, of the lowest vibrational state are 
in reasonable agreement with previous studies. In particular, we confirm a general trend of 
the considerable increase of vibrational frequency viz. growth of B indicated for the first 
time by Larsen ||. The energy dependence on the magnetic field is much more pronounced 
for the lowest rotational state - it grows much faster than the vibrational one with magnetic 
field increase. It implies that the if^ in parallel configuration becomes more stable for larger 
magnetic fields. From a quantitative point of view the results obtained by different authors 
are not in good agreement. It is worth mentioning that our results agree for large magnetic 
fields > 10a. u. with results by Le Guillou et al. |I3| , obtained in the framework of the so 
called 'improved static approximation', but deviate drastically at B — la.u., being quite 



close to the results of Larsen H] and Wille |J11||. As for the energy of the lowest rotational 



state our results are in good agreement with those obtained by other authors (see Table V). 
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TABLE V: Energies of the lowest vibrational (E v n,) and rotational (E ro t) electronic states asso- 
ciated with l g state at 6 = 0°. The indexes in Le Guillou et al |l3| correspond to the 'improved 
adiabatic approximation' (a), and to the 'improved static approximation' (b). 



B 


Et (Ry) 


Evib (Ry) 


Erot (Ry) 




10 9 G 


-1.15070 


0.013 


0.0053 


Present 




— 


0.011 


0.0038 


Wiiic (n) 


la.u. 


-0.94991 


0.015 


0.0110 


Present 




— 


— 


0.0086 


Wiiic (n) 




— 


0.014 


0.0091 


Larsen |9j 




— 


0.013 


— 


Lc Guillou et al (a) Q 




— 


0.014 


0.0238 


Le Guillou et al (b) Q 


10 10 G 


1.09044 


0.028 


0.0408 


Present 




— 


0.026 


0.0308 


Wille Jn| 


10 a.u 


5.65024 


0.045 


0.0790 


Present 




— 


0.040 


0.133 


Larsen [bl 




— 


0.039 


— 


Le Guillou et al (a) Q 




— 


0.040 


0.0844 


Le Guillou et al (b) jl|| 


10 11 G 


35.0434 


0.087 


0.2151 


Present 


100 a.u 


89.7096 


0.133 


0.4128 


Present 




— 


0.141 


0.365 


Larsen [bl 




— 


0.13 


— 


Wunner et al J2l| 




— 


0.128 


— 


Le Guillou et al (a) Q 




— 


0.132 


0.410 


Le Guillou et al (b) jl|| 


10 12 G 


408.389 


0.276 


1.0926 


Present 




— 


0.198 


1.0375 


Khersonskij |q| 


1000 a.u 


977.222 


0.402 


1.9273 


Present 




— 


0.38 


1.77 


Larsen [PJ 




— 


0.39 


— 


Wunner et al J2l[ 




— 


0.366 


— 


Lc Guillou et al (a) Jl|] 




— 


0.388 


1.916 


Le Guillou et al (b) jl|| 


10 13 G 


4219.565 


0.717 


4.875 


Present 




— 


0.592 


6.890 


Khersonskij flq] 


4.414 x 10 13 G 


18728.48 


1.249 


12.065 


Present 



We show the electronic distributions 



[x,z 



integrated over y and normalized to one 



for magnetic fields 10 9 , 10 10 , 10 11 , 10 12 G and different orientations in Fig. 6. It was already 
found [^] that there is a change from ionic (two-peak electronic distribution) to covalent 
coupling (single-peak distribution) at 9 = 0. If for B = 10 9 G, all electronic distributions are 
characterized by two peaks, then for 10 12 G all distributions have a sharp single peak. Fig. 6 
demonstrates also how the change of the type of coupling appears for different inclinations. 
It is quite natural that for the perpendicular configuration 9 = 90°, where the equilibrium 
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distance is the smallest, this change appears for smaller magnetic field. 

To complete the study of the l g state we show in Fig. 7 the behavior of the variational 
parameters of @ as a function of the magnetic field strength for the optimal (parallel) 
configuration, 6 = 0. In general, the behavior of the parameters is rather smooth and very 
slow-changing even though the magnetic field changes by several orders of magnitude. In 
our opinion it reflects the level of adequacy (or, in other words, the level of quality) of our 
trial function. In practice, the parameters can be approximated by splain method and then 
can be used to study any magnetic field strength other than presented here. 

B. 1„ state 



In the absence of a magnetic field the l u state |33|] is essentially repulsive, antibonding and 
characterized by shallow minimum in the total energy situated at large internuclear distance 
(see, for example, H,||' '). It is a weakly bound state with respect to dissociation and 



it becomes unbound if nuclear motion is taken into account. There are not many studies of 
this state. Our major finding is that the total energy surface of the system (ppe) in the state 
l u exhibits a well-pronounced minimum for < B < 4.414 x 10 13 G and for all inclinations. 
Similar to the l g state, both total (Et) and binding (E^) energies of the l u state increase as 
the magnetic field grows, while the equilibrium distance shortens. However, the accuracy of 
our calculations do not allow us to make a definitive conclusion about the stability of the 
system with respect to dissociation and nuclear motion effects. 

Variational results for the l u state at 0°, 45° and 90° for magnetic fields B = 0— 4.414x 10 13 
G are shown in Tables VI- VIII. The immediate conclusion is that Et{0°) < -Er(45°) < 
-&r(90°) for all magnetic fields. Hence, as for the l g state, the highest molecular stability of 
the l u state occurs for the parallel configuration. Also, the binding energy growth is maximal 
as a function of magnetic field for the parallel configuration. Therefore, the stability of B.^ in 
the parallel configuration in the l u state even increases as the magnetic field grows, similarly 
to what happens for the l g state. These results suggest the following picture for appearance 
of bound state: for small magnetic fields the minimum in the total energy arises at very large 
internuclear distances [34], then, as the magnetic field grows, the position of the minimum 



moves to smaller and smaller internuclear distances. 

Our results for B > and 9 = 0° give the lowest total energies compared to other 



22 



B = 10 9 G 




# = 0° 




= 45° 




6» = 90° 



B = 10 1U G 




# = 0° 




= 45° 




90° 



B = 10 U G 



3 

2 


* 


3 
2 


A 


3 
2 


0k 


1 




1 




1 


t ~ \ 





z 1 0.4 x 





z 1 0.4 x 





"' z° 1( X4 




= 0° 




# = 45° 




6» = 90 



£ = 10 12 G 



5 









15 

10 


5 






t*'¥k 


5 





r 


6~ 

z 


10.4 U x 








z u j0.4 u x 

6> = 0° = 45° 

FIG. 6: H2 electronic distributions |^(x,z)| 2 (normalized to one) for the l g state for different 
magnetic fields and inclinations. 
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FIG. 7: Variational parameters of the trial function (8) as function of the magnetic field strength 
B for the l g state in parallel configuration, 6 = 0°. The parameter ^3 is fixed being equal to 1 
(see text). The error bars correspond to relative deviation in the variational energy in the region 
AE T = ^< 10~ 5 . 
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TABLE VI: Total Ex, binding Ef, energies and equilibrium distance R eq for the l u state in the 
parallel configuration, 6 = 0°. ^This value is taken from |E0| 



B 



E T (Ry) 



E b (Ry) 



R eq (a.u. 



B = Q 


-1.00010 


1.00010 


12.746 


Present^ 




-1.00012 


1.00012 


12.55 


Peek et al 


10 9 G 


-0.92103 


1.34656 


11.19 


Present 




-0.917134 


— 


10.55 


Peeket al |17] 


1 a.u. 


-0.66271 


1.66271 


9.732 


Present 




-0.66 


1.66 


9.6 


Kappes et al 14] 


10 10 G 


1.63989 


2.61500 


7.180 


Present 




2.1294 


— 


4.18 


Peeket al 1 17] 


10 a.u. 


6.52362 


3.47638 


6.336 


Present 


10 n G 


36.83671 


5.71649 


4.629 


Present 


100 a.u. 


92.42566 


7.57434 


3.976 


Present 


10 12 G 


413.6175 


11.9144 


3.209 


Present 


1000 a.u. 


984.6852 


15.3148 


2.862 


Present 


10 13 G 


4232.554 


22.7648 


2.360 


Present 


4.414 x 10 13 G 


18750.07 


32.9104 


1.794 


Present 



calculations. They are in good agreement with those by Kappes-Schmelcher |T4|] <a ' ) as well 
as by Peek-Katriel |3j]] for B = 0, 10 9 G, although for B = 10 10 G a certain disagreement is 
observed (see Table VI). However, for 9 = 90° our results are in striking contrast with those 
by Wille [JL2]J , where even the optimal configuration is attached to 9 = 90°, contrary to our 
conclusion. For instance, at B = 10 10 G in |I2|] the values Ej, = 2.593 Ry and R eq = 2.284 
a.u. are given, while our results are J5& = 1.9617 Ry and R eq = 5.517 a.u., respectively (see 
Table VIII). Similar, but less drastic disagreement, is the observed with the results in ||14j| w . 
We can guess this disagreement is due to the shallow nature of the minimum, but a real 
explanation of this effect is missing. 
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TABLE VII: Total Et, binding Ef, energies and equilibrium distance R eq for the l u state in the 
configuration 6 = 45°. Optimal value for the gauge parameter £ is shown (see text). 



B 


E T (Ry) 


E b (Ry) 


R eq (a.u.) 


£ 


10 9 G 


-0.870391 


1.295923 


8.053 


0.9308 


1 a.u. 


-0.509041 


1.509041 


6.587 


0.9406 


10 10 G 


2.267998 


1.987321 


4.812 


0.9671 


10 a.u. 


7.692812 


2.307188 


4.196 


0.9808 


10 11 G 


39.71061 


2.84258 


3.538 


0.9935 


100 a.u. 


96.88464 


3.11536 


3.278 


0.9968 


10 12 G 


422.0074 


3.5245 


3.020 


0.9991 


1000 a.u. 


996.3044 


3.6956 


2.894 


0.9996 


10 13 G 


4251.409 


3.9103 


2.790 


0.9999 


4.414 x 10 13 G 


18778.95 


4.0330 


2.746 


0.9999 



The analysis of Tables VI- VIII shows that for 9 > 0° and fixed magnetic field the total 

energy of H^ in the l u state is always larger than the total energy of the hydrogen atom 

25| . It means that the H^-ion in the l u state is unstable towards dissociation to H + p. 



For 9 ~ 0° the present total energies of the H% ion and the most accurate results for the 
hydrogen atom [T25J are comparable within the order of magnitude 10~ 4 — 10 -5 . At the same 
time we estimate that the accuracy of our calculations is of the same order of magnitude. It 
prevents us from making a conclusion about the stability of H^ in the l u state with respect 
to dissociation. The only conclusion can be drawn is that the minimum is very shallow. 

The l u state is much more extended than the l g state: for fixed magnetic field the 
equilibrium distance of the l g state is much smaller than that for the l u state. This picture 
remains the same for any inclination. It is quite striking to see the much lower rate of 
decrease of R eq in the range 5 = — 4.414 x 10 13 G: for the state l u it falls ~ 3 times 
compared to the l g state, where it falls ~ 20 times. 

The behavior of the equilibrium distance R eq of the l u state as a function of inclination 
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TABLE VIII: Total Et, binding Ef, energies and equilibrium distance R eq for the l u state at 
8 = 90°. Optimal value for the gauge parameter £ is shown (see text). 



B 



E T (Ry) E b (Ry) R eq (a.u.) 



s c 



10 9 G 



1 a.u. 



10 10 G 



0.867234 


1.292766 


8.784 


-0.49963 


1.49963 


7.264 


-0.65998 


1.65998 


5.45 



2.29365 



1.96167 



5.517 



10 a.u. 


7.72998 


2.27002 


4.872 


10 11 G 


39.76500 


2.78819 


4.154 


100 a.u. 


96.93497 


3.06503 


3.875 


10 12 G 


422.0834 


3.44848 


3.594 


1000 a.u. 


996.3807 


3.61935 


3.460 


10 13 G 


4251.497 


3.82238 


3.340 


4.414 x 10 13 G 


18779.04 


3.9409 


3.306 



0.9692 Present 

0.9737 Present 

Kappes et al 0(b) 

0.9866 Present 

0.9923 Present 

0.9975 Present 

0.9988 Present 

0.9997 Present 

0.9998 Present 

0.9999 Present 

0.9999 Present 



is quite non-trivial (see Tables VI- VIII). As in the l g state, the H^~-ion in the l u state for 
B < 10 12 G is most extended in the parallel configuration, while being most compact at 
6 ~ 45°. For B ~ 10 12 G the equilibrium distance is almost independent of the inclination, 
while for B > 10 12 G the most compact configuration corresponds to 9 = 0°, in contrast to 
the l g state. 

IV. CONCLUSION 

We carried out accurate non-relativistic calculations in the Bohr-Oppenheimer approx- 
imation for the lowest states of the H^~ molecular ion of even parity, l g , and odd parity, 
l u , for different orientations of the magnetic field with respect to the molecular axis. We 
studied constant magnetic fields ranging up to B = 4.414 x 10 13 G, where a non-relativistic 
consideration is valid. 

For all studied magnetic fields and orientations a well-pronounced minimum in the total 
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energy surface for both l g and l u states was found. This makes manifest the existence of 
i?2~ in both states for magnetic fields B = — 4.414 x 10 13 G. The smallest total energy was 
always found to correspond to the parallel configuration 9 = 0°. The total energy increased 
while the binding energy decreased steadily as the inclination angle grew from 0° to 90° for 
both states. The rate of the total energy increase as well as the binding energy decrease was 
seen to be always maximal for the parallel configuration for both states. 

However, the equilibrium distance exhibited quite non-trivial behavior as a function of 
the orientation angle 9. In the case of the l g state, the shortest equilibrium distance always 
corresponded to the perpendicular configuration for magnetic fields B < 10 12 G, whereas 
for B > 10 12 it occurred for some angle 9 < 90°. On the contrary for the l u state the 
shortest equilibrium distance always corresponded to orientations 9 < 90° for all magnetic 
fields considered. In particular, for B > 10 12 G, it begins to correspond to the parallel con- 
figuration. As for the largest equilibrium distances - they were found to always correspond 
to the parallel configuration for both the l g state for all studied magnetic fields and for the 
l u state for B < 10 12 . However, for larger magnetic fields B > 10 12 the largest equilibrium 
distance for the l u state was seen to correspond to the perpendicular configuration (!). 

Confirming the previous qualitative observations made by Larsen and Khersonskij |7j 
for the l g state we demonstrated that the H^ ion in the lowest energy state can dissociate 
to H + p for a certain range of orientations starting from magnetic fields B> 1.8 x 10 11 G. 
As the magnetic field increases the region where dissociation is allowed was seen to steadily 
broaden, reaching 25° < 9 < 90° for B = 4.414 x 10 13 G. 

The electronic distributions were found to be qualitatively different for weak and large 
magnetic fields. In the domain B < 10 10 G the electronic distribution for any inclination 
peaks near the position of the protons. On the contrary for B > 10 11 G the electronic 
distribution is always peaked near the midpoint between the protons for any inclination. It 
implies physically different structure of the ground state - for weak fields the ground state 
can be modeled as a combination of hydrogen atom and proton while for strong fields such 
modeling is irrelevant. 

Combining all the above-mentioned observations we conclude that for magnetic fields of 
the order of magnitude B ~ 10 11 G some qualitative changes in the behavior of the H^ 
take place. The behavior of the variational parameters also favors this conclusion. It looks 
like as appearance of a new scale in the problem. It might be interpreted as a signal of a 

28 



transition to the domain of developed quantum chaos 
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